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ABSTRACT 


The thesia ia about algnal analyala uaing wavelet tranaform. 
Uavelet theory, which haa gained a lot of aomentuin haa been 
inveat igat ed . Varioua orthogonal baaia function which make up the 
waveleta have been analyaed. The corr eaponding waveleta have been 
found out. The filtera baaed on the waveleta/orthogonal baaia 
functiona have been deaigned. A computer program haa been 
developed to do the decompoaition and reconatruct ion of aignala 
baaed on an iterative algorithm. The aignala which have been 
taken up for analyaia are compoaite exponentiala and ainuaoidal. 
The aame aignala with cloaely apaced burata and noiae have alao 
been taken up for inveat igat ion. A aignala have alao been 
decompoaed and reconatruct ed by uaing another acheme called 
Laplacian Pyramid Scheme. 

A aeparation of exponentiala , an age old problem, haa alao 
been taken up for inveatigatlon in thia context. The purpoae of 
thia inveatigation ia to aee whether by during wavelet 
decompoaition and reconatruct ion, the exponentiala which 
conatitute the original aignal, aurface or not 
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CHAPTER - 1 


INTRODUCTION 

Uavelet theory, which has of late has eenerated a lot of 
interest, provides a unified framework for a number of techniques 
which had been developed independently for various sifinal 
processinft applications. Nultiresolut ion sianal processinq used 
in computer vision, sub-band coding of speech and image 
compression and wavelet series expansions developed in applied 
mathematics, are all different views of single theory. 

The wavelet theory, which covers quite a large area, treats 
both the continuous and discrete-time cases. It provides very 
general techniques that can be applied to many tasks in signal 
processing and hence, has numerous potential applications. In 
particular, the Uavelet Transform CWT) is of interest for the 
analysis of non-s tat ionary signals, because it provides an 
alternative to the classical Short-Time Fourier Transform CSTFT) 
or Gabor transform. The basic difference between UT and STFT is 
that in contrast to the STFT, which uses a: s4«k«^ u»m a single 
analysis window, the UT uses short windows at high frequencies 
and long windows at low frequencies. 

The Uavelet Transform can be seen as a signal decomposition 
onto a set of basis functions. In fact, basis functions called 
wavelets always underlie the wavelet analysis. They are obtained 
from a single prototype wavelet by dilations and contractions 
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(scalings) as well as shifts. The prototype wavelet can be 
thought of as a bandpass filter. 

In a Uavelet Transform, the notion of scale is introduced as 
an alternative to frequency, leading to a so-called time-scale 
representation. In other words, a signal is mapped into a time- 
scale plane in contrast to the time-frequency plane used in STFT. 

There are several types of wavelet transforms. The choice of 
a transform depends on the type of application. For example, for 
a continuous input signal, the time and scale parameters can be 
continuous, leading to the continuous wavelet transform. They may 
be discrete also, giving rise to a Uavelet Series Expansion. Or 
the wavelet transform can be defined for discrete-time signals, 
which is known as Discrete Uavelet Transform (DUT). In the latter 
case, it uses Hultirate Signal Processing techniques. In this 
context, one can draw an analogy with the Fourier Transform 
(continuous), Fourier Series and Discrete Fourier Transform. 

Uavelet theory has been developed as a unifying framework 
only recently, although similar ideas and constructions took 
place as early as the beginning of the century. The idea of 
looking at a signal at various scales and analysing it with 
various resolutions has in fact emerged independently in many 
different fields. In mid-eighties, researchers such as Morlet, 
Grossmann, and Neyer built stropg mathematical foundations around 


the subject. 
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The attention of the sienal processine community was soon 
caufiht when Daubechies and Hallats, in addition to their 
contribution to the theory of wavelets, established connections 
to the discrete siftnal processins results. 

In the paper presented, various bases of wavelets have been 
invest ifiat ed . The frequency response of the filters correspondina 
to the bases have been explored. In addition, a comparison 
between Uavelet Transform and Short-Time Fourier Transform has 
been drawn. A study has been carried out on the decomposition and 
reconstruction of composite sianals such as exponentials with 
varied time constants and sinusoidals with varied frequencies, 
with the help of the wavelet transform and Laplacian Pyramid 
Scheme. The reconstruction of the above mentioned sianals with 
bursts and noise have also been examined. The separation of 
exponentials, a problem which researchers have been tryine to 
solve usinft various methods, is of sreat significance and 
importance in some fields such as physics. The problem has been 
looked into from the Wavelet Transform point of view. 

In the thesis. Chapter 2 deals with the multiresolution 
transform and its implementation, the wavelet representation, 
pyramidal algorithm, comparison between wavelet trasnform and 
short time fourier transform, and Laplacian pyramid alsorithm. 

In Chapter 3, we have described the results of various 
<l 0 si|[^ns and operations. The chapter also deals with the design of 
wavelets and filters correspondina to different orthoaonal basis 



4 


functions. The results of reconstruction of exponential and 
sinusoidal sifinals have been described. In this chapter, we also 
discuss the problem of separation of exponentials with regards 
to wavelet transform. 

Conclusion has been drawn in Chapter 4. 



CHAPTER - 2 


WAVELET THEORY 


2.1 MOTATIOM 

a) 2 A R denote the set of integera and real numbera reapecti- 

2 

vely. L (R) denotea the vector apace of measurable, square 
integrable one-dlmenalonal functions f(x). 

b) For f(x) « L^(R) and «(x) « L^(R), the inner product of £(x) 
with g(x) is given by 

to 

< g(u), fCu) ^ “ J g(u)-f(«) dw 

-to 


c) The norm of f(x) is L (R) ia given by 

to 


li f ir 


= I |fCu)| 


du 


d) The convolution of two functions f(x) « L (R) and g(x) 
L^(R) ia denoted by 


f * g(x) = * «(«)] (X) 

a J f(u).g(X-U) du 


e) The Fourier Transform of f(x) « L (R) is written fCw) and is 
defined by 

to 

£(w) ® J f(x) exp(-wx) dx 
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(£) I (Z) is the vector space of square sustmable sequences 

2.2 MULTI RESOLUTl CM TRANSFOAM 

2.2.1 MultlrAsoluilon Approximation of' L^CK> [2] 

Let A^J he the operator which approximates a signal at 

resolution 2^. Ue suppose that our original signal f(x) is 

2 

measurable and has a finite energy : f(x) « L (R) . Here, we 
characterize A^j from the intuitive properties that one would 
expect from such an approximation operator. Ue state each 
property in words, and then give the equivalent mathematical 
formulation. 

(i) A^j is a linear operator. If is the approxioiation 

of some function f(x) at the resolution 2^; then A 2 jfCx) i« 

not modified if we approximate it again at the resolution 

2-^. This principle shows that A 2 j i« a projection operator 

2 

on a particular vector space V 2 j <= L (R)- The vector space 
V.j can be interpreted as the set of all possible 

m 

3 2 

approximations at the resolution 2-' of functions in L CR) - 

(ii) Among all the approximated functions at the resolution 2 ^, 
A„jf(x) is the function which is the most similar to fCx). 

<4 


V g(x) « V2j, llg(x)-f(x)ll a llA2j£(x) - £(x)|| 


( 2 . 1 ) 
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Hence, the operator A 2 J ia an orthogonal projection on 
the vector apace ^ 2 ^' 

(ili) The approximation o£ a aignal at a reaolution 2 ^*^ containa 
all the neceaaary Information to compute the aame aignal at 
a amaller reaolution 2 ^. Thia ia a cauaality property. 
Since A^j ia a projection operator on V^j, thia principle 
ia equivalent to 

Vj « Z, V^j c V2j+1 C2.2) 

(iv) An approximation operation ia aimilar at all reaolutiona. 
The apacea of approximated functiona ahould thua be derived 
from one another by acaling each approximated function by 
the ratio of their reaolution valuea. 

Vj « Z, f(x) € V^j <=> f(2x) « V^j+l (2.3) 

(v) The approximation A^jfCx) of a aignal £(x) can be 

characterized by 2-^ aamplea per length unit. Uhen f(x) ia 
tranalated by a length proportional to 2 A 2 jf(x) ia 

tranalated by the aame amount and ia characterized by the 
aame aamplea which have been tranalated. Aa a conaequence 
of (2.3), it la aufflclent to expreaa the principle (V) for 
the reaolution j^o. The mathematical tranalationa conaiat 
of the following 

. Dlscrmta characterization s 

2 

There exiata an iaomorphiam I from onto I (Z) (2.4) 
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. Translation of the approximation i 


Vk e 2, A^£j^(x) = A^fCx-k). where £j^(x) = £(x-k) (2.5) 

Tranalation o£ the aamplea : 

<=> I ^ <2-« 


(vi) Uhen computing an approximation o£ £(x) at resolution 2^, 
aome information about £(x) ia loat. However, aa the 
resolution increase to +<n, the approximated signal should 
converge to the original signal. Conversely aa the 
resolution decreases to aero, the approximated signal 
contains leas and less Information and converges to aero. 


Since the approximated signal at a resolution 2^ is 
equal to the orthogonal projection on a space ^2-^’ this 
principle can be written 

Lim V.j = V-i is dense in L^(R) (2.7) 
j -► +o» ^ j=-«> 


and 


Lim 

j -► -<» 


^2^ 




^ 2 ^ 


Ue call any set of vector spaces 


{“} 


( 2 . 8 ) 


which satisfies 


the properties (2.2)-(2.d) a multiresolution approximation of 

L^(R). The associated set of operators A^j satisfying (i) and 

2 

(vi) give the approximation of any L (R) function at a resolution 
zK Ue now give a simple example of a multiresolution 


approximation of L (R) 
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2 

ixampla i Let be the vector space of all functions of L (R) 
rhich are constant on each interval ]K,K+1[ for any k«Z. Equation 
;2.3) implies that V^j is the vector space of all the functions 
>f L^(R) which are constant on each interval ]k2”^, (k;+l)2~-^[, 

for any k « Z. The condition (2.2) is easily verified. He can 
lefine an isophormism I which satisfies properties (2.4), (2.5) 

md (2.6) by association with any function f(x) c the sequence 
C«k^k«Z that equals the value of f(x) on the interval ]k, 

k:+l[. Ue know that vector space of piecewise constant functions 

2 QD 

is dense in L (R). Hence we can derive that . u V.,j is dense in 

j=-® a-* 

L^(R). It is clear that sequence of vector 

spaces fv_jl is a multiresolution approximation of L^(R). 
^ Jj«Z 

Unfortunately, the functions of these vector spaces are neither 
smooth nor continuous, making this multiresolution approxi- 
mation rather inconvenient. 

Ue saw that the approximation operator A^j is an orthogonal 
projection on the vector space In order to numerically 

characterize this operator, we must find an orthogonal basis of 
V^j. The following theorem shows that such an orthogonal basis 
can be defined by dilating and translating a unique function 
^(x). 

ThaoTMi 1 I Let fv.,j| be a multiresolution approximation of 
“I 2 J j«Z 2 

L (R) . There exists a unique function <^(x) e L (R) , called the 
scaling function, such that if we set “ 2V(2^x) for j«Z, 

(the dialation of ^(x) by 2^), then 
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is an orthogonal baais oi V . (2.9) 

The theorem ahowa that we can build an orthonormal baala o£ 
any V^j by dilating a function <^(x) with a coefficient 2^ and 
tranalating the reaulting function with a grid whoae interval ia 
proportional to 2 . The function ^jjjCx) are normalized with 

reapect to L^(R) norm. The coefficient /Tf appeara in the baaia 
aet in order to normalize the functiona in the L^(R) norm. For a 
given mult ireaolution approximation (V j) . , there exiata a 
unique acaling function 4^(x), which aatiafiea (2.9). However, for 
different mult ireaolution approximationa , the acaling functiona 
are different. 

The orthogonal projection on ^ 2 ^ computed by 

decompoaing the aignal f(x) on the orthogonal baaia given by 
Theorem 1. Specifically, 

V f(x) « L^(R), 

to 

A 2 j£(x) = 2”-^ ^ < f(u), ^ 2 ^ (u-2"jn) > 4,^1 (x-2"jn) (2.10) 

The approximation of the aignal f(x) at the reaolution 2^, 
A 2 J£(x), ia thua characterized by the aet of inner producta which 
we denote by 

A^j £ . [< f(u). ♦2 JCu-2'Ja) >]^ 


( 2 . 11 ) 
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la called a diacrete approxlttation of f(x) at the reaolutioft 
2^. Each Inner product can alao be interpreted aa a convolution 
product evaluated at a point 2~^n 


w 

<£(u) > * J £(u) .^2 j(^~2"^n) du 

= (f(u) * ^ 2 jC-u)) (2"-*n) 
Hence, we can rewrite A^jf: 

A^Jf = [Cf(u) « [2'^")]n.z 


( 2 . 12 ) 


Since ^(x) ia a lowpaaa filter, thia diacrete aifinal can be 
interpreted aa a lowpaaa filtering of f(x) followed by a uniform 
aampling at the rate 2^. In an approximation operation, when 
removing the detaila of f(x) amaller than 2 we auppreaa the 

higheat frequenciea of thia function. The acaling function ^(x) 
forma a very particular lowpaaa filter aince the family of 

^2j(x-2‘-’n) 

2.2.2 Implementation of a Multiresolutlon Transform [21 


irmZ 


is an orthonormal family. 


functions 




In practice, a physical measuring device can only 
measure signal at a finite resolution. For normalization 
purposes, we suppose that thia resolution is equal to 1. Let 
be the discrete approximation at the resolution 1 that ia 
measured. The causality principle says that from A^f we can 
compute all the discrete approximations A^j for j<0. In this 
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section, a simple Iterative algorithm Is described for 
calculating these discrete approximations. 

Let * nultlresolutlon approximation and <jfr(x) be 

the corresponding scaling function. The family of functions 

/ mmm j mm " l 

2 4 ^ 2 ^ orthonormal basis of V^j + l. Ue 

know that for any n « Z, the function 4 ^ 2 jCx-2~-^n) Is a member of 
^2'^ which Is Included In V 2 j + 1- It can thus be expanded In this 
orthonormal basis of V 2 j +1 : 

^23(x“ 2 ^n) = 

00 

2"-*“^ ^ < ^23(«-2"^n) ^2i + lCu-2" ^'^k) > ^^’} + lCx-2~ 

K=-® 

(2.13) 

By changing variables In the Inner product Integral, one can 
show that 

2"-^“^ < ^2iC’*-2”-^n)- 4fr2j+l<«-2"-^"^k) > 

= < ^(u-(k- 2 n)) > 

Uhen computing the Inner products of f(x) with both sides of 
(2.13), we obtain 

<n 

< f(u), ^ 2 jCu- 2 "-’n) > » ^ 4 fr(tt-(k- 2 n)) > 

k*=-oo 

< f(u), ^23+l(«-2"-*"^k) > (2.14) 
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Let H be the discrete filter whose impulse response is given 
by 

Vn « Z . h(n) = < ^(u-n) > (2.15) 


Let H be the mirror filter with impulse response h(n) * h(“n). By 
inserting (2.15) in the previous equation, we have 


CO 


< f(u), (u-2"^n) > = ^ h (2n-k) <£(u), + *k) 


k=-<B 


(2.16) 


Equation (2.16) shows that ^2^^ computed by the 

convolution of A^j+lf and H and keeping every other sample of the 

di 

output. All the discrete approximations A^jf* £or j < 0; can thus 
be computed from A^f by repeating this process. This operation is 
called a pyramid transform. The algorithm is illustrated by a 
block diagram later in the chapter (Fig. 2.1). 


In practice, the measuring device gives only a finite number 

of samples : A?f * (« ) ^ m- Each discrete signal A^jf(j<0) 

1 nlsnSN 4 

haa 2 *^N aamplaa* In order to avoid border problema, when 
computing the diacrete approximat iona we auppoae that the 

original signal A^f is symmetric with respect to n»0 and n»N. 

{ oi_^ if -N < n < 0 

-2N-n 0 < n < N 

If the impulse response of the filter H is even (H=H) , each 
diacrete approximation A^jf will also be symmetric with respect 
to n=0 and n=2 
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Theorem 1 ahowa that a multireaolut ion approximation ytiZ 

is completely characterized by the scaling function, ^(x) . A 

2 

scaling function can be defined as a function ^(x) « L (R) such 
that, for all j « Z, (/^ orthonormal 

family, and if V^j is the vector space generated by this family 
of functions, then approximation of 

L^(R). 


Theorem a 1 Let 4^(x) be a scaling function and let H be a 
discrete filter with impulse response h(n) = 0(u-n)>. 
Let H(w) be the Fourier Transform defined by 


H(6>) = ^ h(n) e“^“" 

n=-«> 

H(w) satisfies the following two properties : 
H(0) = 1 

and 

|H(6>)|^ + |H(w + n)|^ 

Conversely, |H(6>)| s»* 0 for <i> « [0, n/2] 

The function defined by 

p 

4 ^ («) “ n H(2 <>>) 

p=l 

is the Fourier transform of a scaling function. 


(2.17) 


C2-17a) 


(2.17b) 

(2.17c) 


(2.18) 


The filters that satisfy property (2.17b) are called 
conjugate filters. Given a conjugate filter H which satisfies 
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|[2.17a) - (2.17c), we can then compute the Fourier transform of 
the corresponding scaling function with (2.16). One can easily 
show that the corresponding function H(6>) satisfies 

H(<iJ) = exp(-j6>) co8(<o/2) 

a. 3 THE WAVELET REPRESENTATION 

In this section, multiresolut ion representation based on the 

differences of the information available at two successive 
i i + 1 

resolutions 2'' and 2*' is developed. It is shown that such a 
representation can be computed by decomposing the signal using a 
wavelet orthonormal basis. 

3.3.1 Thm Detmll Signal [21 

Here, we explain how to extract the difference of 

information called the detail signal at the resolution 2-^. The 

HI i 

approximation at the resolution 2'' and 2 of a signal are 
respectively equal to its orthogonal projection on V^J+l and ^ 2 ). 
By applying the projection theorem, it can be easily shown that 
the detail signal at the resolution 2-^ is given by the orthogonal 
projection of the original signal on the orthogonal complement of 
V, j in V„j+1. Let 0,j be this orthogonal complement l.e. 

02 j ia orthogonal to V 2 j, 

“2J • ''2J * ''2^** • 

To compute the orthogonal projection of a function f(x) 
on 02 j, it is required to find an orthonormal basis of 02 j. 
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Theorem 3 ehowa that auch a haaia can be built by acaling and 
translating a function yCx). 

Theorem 3 i Let * multiresolution vector space 

sequence, ^(x) the scaling function, and H the corresponding 
conjugate filter. Le^ V>(x) be a function whose Fourier transform 
is given by 

A M 

V(6>) = G(u/2) .^(6i/2) 

with G(<i>) «s e"’*^ Hfw+n) (2.19) 

Let " 2'^y'(2'^x) denote the dilation of vCx) by . Then 

“ orthonor.al b.ai. of O^j and 
(.^2 ^ ^^jCx-2 ■^*1))^-^ j)eZ orthonomal basis of L^(R). y'(x) 

is called an orthonomal wavelet. An orthonormal basis of O^j can 
thus be computed by scaling the wavelet yCx) with a coefficient 
2-^ and translating it on a grid whose interval is proportional to 
2 "^. 


The decomposition of a signal in an orthonormal wavelet 
basis gives an intermediate representation between Fourier and 
spatial representations. Due to this double localization in the 
Fourier and Spatial domains, it is possible to characterize the 
local regularity of a function f(x) based on the coefficients in 
a wavelet orthonomal basis expansion. 


° 2 ->- 


Let be the orthonormal projection on the 

®2 

As a consequence of Theorem 3, this operator 


vector space 
can now be 


written 
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OOr 




■* ^ < £(u), >.W 2 J(*- 2 “^n) (2.20) 




Pq j^(x) yielda to the detail aignal of f(x) at the resolution 2^. 
2 

It is characterized by the set of inner products 


>J^2 (2.21) 

D^jf iB called the Discrete Detail Signal at the resolution 2'^. 
It contains the difference of information between A^j+lf and 
. It can be shown that each of these inner products is equal 
to the convolution of f(x) with evaluated at 2 "^n 

< f(tt), V 2 J(«- 2 “^n) > = [f(u) * V 

Equations (2.21) and (2.22) show that the discrete detail signal 
at the resolution 2^ is equal to a unifomi sampling of (f(u) * 

V»_j(-u)(x)) at the rate 2^ 

tk 


j(-u)j [2“^»j (2.22) 


°2^' 


» (f(u) » 


¥'2jC-'i) 


(2 




ricZ 


The wavelet V'(3c) can be viewed as a bandpass filter whose 
frequency bands are approximately equal to [-2fr, -fi ] U (n , 2ft]. 

Hence, the detail signal D^jf describes f(x) in the frequency 
bands |^-2”-^'*’^ft, -2~^fij u ^2“-’n, 2 ■ 

It can be proved that by induction that for any J > 0, the 
original discrete signal A^f measured at the resolution 1 is 
represented by 

[aJ-ji. s J s -1 ) 


(2.23) 
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Thia aet of dlacrete aignala ia called an orthogonal wavelet 
repreaentat ion , and conaiata of the reference aignal at a coarae 
reaolution A^'jf and the detail aignala at the reaolutiona 2"^ for 
-J S j i -1. It can be interpreted as a deconpoait ion of the 
original aignal in an wavelet baaia. 

2.3.13 Implementatrlon of an Orthogonal Wavelet. Representation 131 

In thia aection, we deacribe a pyranidal algorithm to 
compute the wavelet repreaentat ion. 

For any rieZ, the function "^n) is a member of 

O^j In the same manner as (2.13), this function can be 

expanded in an orthonormal basis of V^j+l. 

CO 

= 2"-^"^ ^ (u-2"^n), * 2 ^*^ (u-2"J"^k) > 

k»-<B 

. 4 ^ 2 ^*^ (x-2“-^"^k) (2.24) 

Aa we did in (2.14), by changing the variables in the inner 
integral we can prove that 

2 ^ ^ < ^^iin-Z ^n), 4t^’i+lin-Z ^k) > 

= < ^(«“(lt-2n)) > (2.25) 

Hence, by computing the inner product of f(x) with the functions 
of both sides of (2.24), we obtain 
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< f(u), ^n) > = 

< V' 2 ~ 1 C«)- <^(u-(k-2n)) 

k=-<» 



>.<£Cu),^ 2 j+lCu-2“ j“^k)> 


(2.26) 


-fet G b* the difitct’ete filter with iatpulee response. 
«(n) = < ^(u-n) > 


(2.27) 


3 be the mirror filter with impulse response e(n) = e(-n). 

Inserting (2.27) to (2.26) yields 


< f(u), > = ^ 2 ^^ 


00 


^ g (2n-k) <f(u), ^,j+l (u-2"^"^k)> 


k=-<n 


(2.28) 


Equation (2.28) shows that we can compute the detail signal ^ 2 ^^ 

by convolving A^j+lf with the filter G and retaining every other 

sample of the output. The orthogonal wavelet representation of a 

discrete signal A^f can therefore be computed by successively 
d d 

decomposing A 2 j + lf into A 2 jf and D 2 jf for -J aJ j a: -1. 

s 

Equation (2.19) of Theorem 3 implies that the impulse 
response of the filter G is related to the impulse response of 
the filter H by 

g(n) = (-1)^"“ h(l-n) (2.29) 


G is the mirror filter of H, and is a highpass filter. In signal 
processing, G and H are quadrature mirror filters. Equation 
(2.28) can be Interpreted as a highpass filtering of the discrete 
signal A^j+li- 
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If the original signal has N samples, then the discrete 

signals D^jf and A^jf have 2'^N samples each. Thus, the wavelet 

c3> 

representation (A 2 “jft (D^jf) - J ^ j i -1) has the same total 

number of samples as the original approximated signal A^f. This 

occurs because the representation is orthogonal. The energy of 

the samples of D^jf gives a measure of the irregularity of the 

1 + 1 

signal at the resolution 2-' . Uhenever A^j(x) and A^j + lf(x) were 

significantly different, the signal detail has a high amplitude. 
In Fig. 2.1, the implementation of wavelet decomposition of is 
shown . 


Z. 3. 3 Signal Steconstructlon tram an Orthogonal Wavelet 
Representation 


Ue have seen that the wavelet representation is 
complete. Ue now show that the original discrete signal can also 
be reconstructed with a pyramid transform. Since ^ 2 ^ the 

orthogonal complement of V^j in V^j+l, 

orthogonal basis of V 2 j + 1. For any n > 0, 

— 5-1 

the function ^-j + l(x - 2 n) can thus be decomposed in this 

*1 

basis 


^^ 3*1 [x - 


to 


2 “^ ^ (u-2"^k), 




00 


2 “^ ^ < v'2J('»"2“^k), ^ 2 ^*^ (u-2“^"^n) > .^2 J Cx-2” ^k) 


k=-® 


(2.30) 



u a 



Fig. 21 : D*conpo«ition o£ a discrete approximation A^j+l into an 
approximation at a coarser resolution A^j £ and the 
signal detail £- By repeating in cascade this 

algorithm £or -1 i j ^ -J; we compute the wavelet 

representation o£ a signal A^j £ on J resolution 
levels. 
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By computine the inner product o£ each aide o£ equation (2.30) 
with the £unction £(x) we have 

< f(u), * 2 ^*^ (u- 2 "-*'^) > 

= 2 ^ ^ < ^ 2 j(u -2 ^k), > . <f (u) , 4^2 J ('*" 2 " ^k) > 

k^-oo 

00 

+ 2 "^ ^ < V'gj (u- 2 "-’k), ^2 j + lC'»“ 2 ~^”^n)> <f(u).V 2 J(«- 2 "-^k)> 

k=-o 6 

(2.31) 

Inaerting (2.14) and (2.25) in thia expreaaion and uaing the 
£iltera H and G, reapectively , de£ined by (2.15) and (2.27) 
yielda 

<(f(u), ^ 2 ^*^ ■ 2 "^"^n)> 

= 2 ^ h(n- 2 k).<f(u), (u- 2 "jk)> 

k=-oft 

«o 

- I - n-2k) <£(u), (u-2~^k)> (2.32) 

k=-oe 

Thia equation ahowa that A^j+l can be reconatructed by 
putting zeroa between each aample o£ and and convolving 

the reaulting aignala with the £iltera H and G, reapectively . The 
block diagram in Fig- 2.2 illuatratea thia algorithm. The 
original diacrete aignal A^£ at the resolution 1 is reconatructed 



Ili2> & 




Fie- Reconatructloft o£ a dlacrete approximation A^J+i £ £rom 

41 ^ 

an approximation at a coaraer reaolution A^j £ and the 
aignal detail D 2 j £- By repeating in caacade thia 
algorithm £or -J S j S -1, we are recoatructed £rom ita 
wavelet repreaentation. 
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by repeating this procedure for -J < j < 0. From the discrete 
approximation A^f, we can recover the continuous approximation 
Aj^£(x) with equation (10). 

2.4 MULTI RESOLUTION PYRAMID [41 


Given an original sequence x(n), n<eZ, we derive a lower 
resolution signal by lowpass filtering with a half band lowpass 
filter having impulse response g(n). Following Nyquist rule, we 
can subsample by two (drop every other sample), thus describing 
the scale in the analysis. This results is a signal y(n) given by 


y(n) 


■ I 

Va-< 


«(k) 


x(2n-k) 


The resolution change is obtained by the lowpass filter 
(loss of high frequency detail). The scale change is due to the 
subsampling by two, since a shift by two in original signal x(n) 
results in a shift by one in y(n). 


Now, based on this lowpass and subsampled version of x(n), 
we try to find an approximation, a(n) , to the original. This is 
done by first upsampling y(n) by two (that is, inserting a zero 
between every sample) since we need a signal at the original 
scale for comparison. 

y'(2n) » y(ft) . y'(2n+l) = 0 

Then, y’(n) is interpolated with a filter with impulse 
response g’(n) to obtain the approximation a(n) 
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QO 


a(n) 



a’(k) y*(n-k) 


(2,33) 


k=“0o 

Note that 1£ a(n) and s’(n) were perfect halfband filters 
(having a frequency passband equal to 1 over the nornalised 
frequency range -r/2 to n/2 and equal to 0 elsewhere), then the 
Fourier transform of a(n) would be equal to the Fourier transform 
of x(n) over the frequency range (-n/2, n/2) while being equal to 
aero elsewhere. That is, a(n) would be perfect halfband lowpass 
approximation to x(n) . 


Of course, in general, a(n) is not going to be equal to 
x(n) . Therefore, we compute the difference between a(n) (our 
approximation based on y(n)) and x(n). 


d(n) * x(n) - a(n) 

It is obvious that x(n) can be reconstructed by adding d(n) 
and a(n) , and the whole process is shown in Fig. 2.3. 

In case of a perfect halfband lowpass filter, it is clear 
that d(n) contains exactly the frequencies above n/2 of x(n) , and 
thus, d(n) can be subsampled by two as well without loss of 
information. 

The scheme can be iterated on y(n) creating a hierarchy of 
lower resolution signals at lower scales. Because of that 
hierarchy and the fact that signals become shorter and shorter, 
such schemes are called signal pyramids. 




Pyramid achema. Derivation of lowpaaa, aubaampled 
approximation y(n), from which an approximation a(n) to 
x(n) is derived by upsampling and interpolation. Then, 
the difference between the approximation a(n) and the 
original x(n) la computed as d(n). Perfect reconstruction 
is simply obtained by adding a(n) back. 
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.5 SUB-BAND CODING SCHMES [4] 

The above syateia (pyramid scheme) creates a redundant set o£ 
lamples, since a signal with sampling rate is mapped to two 
tignals d(n) and y(n) with sampling rates £ and £ /2 
i'espect ively . More precisely, one stage o£ a pyramid 
lecomposit ion leads to both a hal£ rate low resolution signal and 
i £ull rate di££erence signal, resulting in an increase in the 
lumber o£ samples by 50%. 

Now a different scheme instead, where no such redundancy 
occurs is looked into. It is the so-called sub-band coding scheme 
first popularized in speech compression. The lowpass, subsampled 
approximation is obtained exactly as explained above, but, 
instead of a difference signal, we compute the "added detail” as 
a highpass filtered version of x(n) (using a filter with impulse 
response h(n)) followed by subsampling by two. Intuitively, it is 
clear that the added detail to the lowpass approximation has to 
be a highpass signal, and it is obvious that if 6(n) is an ideal 
halfband lowpass filter, then an ideal halfband highpass filter 
h(n) will lead to a perfect representation of the original signal 
into two subsampled versions. 

This is exactly one step of a wavelet decomposition using 
sin(x)/x filters; since the original signal is mapped to a 
lowpass approximation (at twice the scale) and an added detail 
signal (also at twice the scale). In particular using these ideal 
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filters, the discrete version is identical to the continuous 
vavelet transform. 

It is not necessary to use ideal (hence impractical) filters 
and yet x(n) can be recovered from its two filtered and 
subsampled versions which we now call yjj(n) and yj^(n). To do so, 
both are upsampled and filtered by g’(n) and h’(n) respectively, 
and finally added together, as shown in Fig. 2.4. Now unlike the 

A 

Pyramid case, the reconstructed signal (which we now call x(n) 
is not identical to x(n) unless the filters meet some specific 
constraints. Filters that meet these constraints are said to have 
perfect reconstruction property. The easiest case to analyse 
appears when the analysis and synthesis filters are identical and 
when perfect reconstruction is achieved i.e. x(n) = x(n) within a 
possible shift. Then it can be shown that the sub-band analysis/ 
synthesis corresponds to decomposition onto an orthonormal basis, 
followed by a reconstruction which amounts to summing up the 
orthgonal projections. FIR filters are assumed. Then it turns out 
that the highpass and lowpass filters are related by 

h(L-l-n) = (-1)*^ g(n) (2.34) 

where L is the filter length. The modulation by (-1)“ transforms 
the lowpass filters to the highpass one. 

2.6 LAPLACIAN PYRAMIOAL MULTI RESOLUTIOM DECOMPOSITION AND 

RECONSTRUCTION 111 

The approximation of a signal f(x) at a resolution ’r’ is 
defined as an estimate of f(x) derived from 'r' measurements per 





J*. 

KCn 


Ca!> 



Fis- 2.H Suband Codlne acheme. (a) Two 

ona corraapondine to low and tha othar to ** S 

ciaa, ara coaputad. Tha caconatructad aignal ia obtainad 

by ra-lntarpolatine tha approxiaationa and 

The flltana on the left £ota an analyaia niftck 

while on tha right ia a aynthaaia 

diagraa (Filter Bank tree) of the 

Tranafora iaplaaantad with diacrete-tiaa filtara 
aubaaapling by two. 
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unit length. These ffleasurements are computed by uniforinly 
sampling at a rate r the function £(x) smoothed by a lowpass 
filter whose bandwidth Is proportional to r. In order to be 
consistent when the resolution varies, these lowpass filters are 
derived from a unique function d(x) which is dilated by the 
resolution factor r : 6^ - •/r&Crx') . The set of measurements 
Arf = (f * called the discrete approximation of 
f(x) at the resolution r. The discrete approximation of a 
function fCx) at the resolution 2^ is thus given by 


(tt]). 


(2.35) 


In pyramidal multiresolution algorithms, the lowpass filter 
function d(x) is chosen such that its Fourier transform can be 
written 

d(e>) » n U (2.36) 

p=l ^ ^ 

where U(e”^^) is the transfer function of a lowpass discrete 

filter U = fu I Let us suppose that we have already computed 

I nj<^4r 

2 

the discrete approximation of a function f(x) L (R) at the 

resolution 2^"^^ ; A 2 j+ 1 £= discrete 

approximation of f(x) at a resolution 2-^ is calculated by 

filtering A 2 j + lf with the discrete lowpass filter U = ['^njrvcZ 

keeping every other sample of the convolution product. 

Let A = I X I „ be such that 
^ ^ J n4i4 

A = A2j+lf * U 


(2.37) 
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th.n AjJ f =. 


(2.38) 


The approximation o£ thia signal at any resolution 2 ~^ , J>0, 
can be computed by iterating on (2.37) and (2.38) and j varying 
between 0 and J+1. This pyramidal algorithm is illustrated in 

OStjSt-J 

called a Gaussian Pyramid. 


Fig. 2.5(a). The sot of discrete approximations 


Now the algorithm by Burt and Crowley is described in order 
to extract the details of f(x) which appear in A^j+lf but not in 
A^jF- The discrete approximation A^j'^'lf has twice as many samples 
as A^jF, so we first expand A^Jf by a factor of two. This is done 
by interpolation. A zero is put between each sample of ^ 2 ^^ 
filter the resulting signal with a lowpass filter. In this 
algorithm, the lowpass filter is the filter U defined previously. 
Let A^jf be the expanded discrete signal. The details at the 

resolution 2-' are then computed by subtracting A 2 jf from A 2 j+ 1 £- 

DaJf = A2j+1£ - A^jf (2.39) 

This algorithm decomposes a discrete approximation A^^f at the 
resolution of 1 into an approximation ^ 2 ”^ ^ 

resolution 2 ’^ and the successive detail signals j^-J’ 

the signal Aj^f has N non-zero samples, each detail signal D 2 jf 
has 2‘^^^N samples; whereas the coarse signal A 2 -j£ has 2 
samples. Hence, the total number of samples of this representa- 
tion is approximately 2N. The signals {V«- ‘''2J«0<JS-J 

regrouped in a data structure called a Laplacian Pyramid. 


are 



♦ 


Fig. 



4 * ‘ 


* 


ca:> 





cw 



k«ap OMr Banipl* out at two 



I 


put » z*r-o botwoon ooch Bampl* 



convolution with tho discrete filter U 


2*5 (a) Decottposltlon of A 2 j + 1 f into A 2 J f and D 2 J f when 

conputina a Laplacian pyraotld (b) Reconstruction of 

A-J+1 f from A-j f and D,j f when reconstructing the 
2 4 

original signal fros a Laplacian Pyramid Scheme. 
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The original aignal can easily be reconstructed froia such a 
decomposition. At each resolution, we compute A2j+1£ by expanding 
the details D 2 j£. By repeating this algorithm when j is varying 
between -J and 0, we reconstruct A^^f. The reconstruction 
algorithm is illustrated by a block diagram in Fig. 2.5Cb) [1]. 

a. 7 COMPARISON BETWEEN SHORT TIME FOURIER TRANSFORM AND WAVELET 
TRANSFORM 

a. 7.1 Short Tims Fourier Transform CSTFT> as Window Fourier 
Transform til {Al 

From the Fourier Transform of a function f(x), we get a 
measure of the irregularities (high frequencies) but this 
information is not spatially localized. Indeed, the Fourier 
transform f(to) is defined through an integral which covers the 
whole spatial domain. It is therefore difficult to find the 
position of the irregularities. In order to localize the 
information provided by the Fourier Transform, Fabor defined a 
new decomposition using a spatial window g(x) in the fourier 
integral. This window is translated along the spatial axis in 
order to cover the whole signal. At a position u and for a 
frequency <■>, the window fourier transform of a function f(x) e 
L^(R) is defined by 

Gf(«,u) * J g(x-u) f(x) dx. (2.40) 
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I£ measurea, locally around the point u, the amplitude o£ 

the sinusoidal wave component of the frequency 6>. In the orifiinal 

Gabor transform, the window function e(x) ia a Gauaaian. It has 

seen been generalized for any type of window function and is 

called a window Fourier Transform. The window function is 

generally a real even function and the energy of its Fourier 

Transform is concentrated in the low frequencies [1]. It can be 

viewed as the impulse response of a lowpass filter. A window 

Fourier Transform can also be interpreted as the inner products 

of the function f(x) with the family of functions 

2 

(g<.>,u(x))(64,u)«R . 


Gf(«,u) = < £(x), g (X) > 

y U 


(2.41) 


In quantum physics, such a family of functions is called a 
family of coherent states. The Fourier Transform g (x) is 

«,U 

given by 


g (<■>) = 

6>_U« 

0 0 


-iUo“ 


g(t4 - Wq) 


(2.42) 


where g(<i>) is the Fourier transform of g(x). A family of coherent 
states thus corresponds to a translation in the spatial domain 
(parameter u) and in the frequency domain (parameter u) of the 
function g(x). This double translation is represented in phase- 
space where one axis corresponds to the spatial parameter u and 
the other to the frequency parameter « (Fig. 2.6). Familiies of 
coherent states have found many applications in quantum physics 
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because they make it possible to analyze simultaneously a 
physical phenomena in both the spatial and frequency domains. 

Let us now describe how a window Fourier transform relates 
to a spatial or a frequency representation. Let be the 

standard deviat ior^iof e(x) 

+0& 

% “ J ^ (2.43) 

Let be the standard deviation of the Fourier transform of e;(x) 

+00 

“ J ^ (2.44) 

“00 

The function g (x) is centered in u^ and has a standard 

“o'"o " 

deviation of in the spatial domain. Its Fourier transform 
u 

£iven by (2.42) is centered in oo^ and has a standard deviation 
By applying the Parseval theorem on (2.41), we get 

+00 

“00 

+0O 

= f f(<^) - (2.45) 

-to 

The first integral shows that in the spatial doma,in, Gf(<«>Q,UQ) 
essentially depends upon the values of f(x) for x «^UQ-o>^,UQ+«»^j . 
The second integral proves that in the frequency domain, 
Gf(<j„,q^) depends upon the values of f(<.>) for 1 . 
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The spatio-frequency domain which ia covered by GfCfo^.UQ) can 
thua be repreeented in the phaae-apace by the reaolution cell 
%] * [“O ■ shown in Fi«. (2.6). The 

surface and shape of the resolution cell is independent from Ug 
and <i>g . The uncertainty principle applied to the function fiCx) 
implies that 

The reaolution cell can therefore not be smaller than 2y2n. The 
uncertainty inequality reaches its upper limit if and only if gCx) 
is a Gaussian. Hence, the resolution in the phase-space is 
maxifflimized when the window function is a Gaussian as in the Gabor 
transform. 


The STFT maps the signal into a two-dimensional function in 
a time-frequency plane (6>,u). It can also be viewed as follows. 
At a given frequency the STFT amounts to filtering the signal 
'at all times’, with a bandpass filter having as impulse response 
the window function modulated to that frequency. Given a window 
function g(x) and its Fourier transform G(<i>), define the 
bandwidth Au of the filter as 


2 |G(6>)|^ d<^ 

® 

J |G(«)P d« 


(2.47) 


where the denominator is the energy of the g(x). Two sinusoids 
will be discriminated if they are more that Au apart. Thus, the 
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resolution in frequency is given by Af . Similarly, the spread in 
time is given by Ax as 


Ax 


Jx^ |g(x:)|^ dx 


(2.48) 


J lg(x)| dx 

where the discrimator is again the energy of g(x). Two pulses in 
time can be discriminated only if they are more than Ax apart. 


Now resolution in time and frequency cannot be arbitrarily 
small, because their product is lower bounded. 

Time-bandwidth product » Ax Au ^ 

This is referred to as the uncertainty principle or 
Heisenberg inequality [4]. It means that one can only trade time 
resolution for frequency resolution, or vice-versa. Gaussian 
windows are therefore often used since they meet the bound with 
equality. 


2.7.1 Wavelet Transform 111 


Horlet defined the wavelet transform by decomposing the 
signal into a family of functions which are the translation and 
dilation of a unique function wC*) ■ The function w(x) ia called a 
wavelet and the corresponding wavelet family is given by 
(-/iV'(«(x-«) ) o2- The wavelet transform of a function f(x) « 
L^(R) is defined by 


00 

Uf(a,u) = J 
*00 


f(x) Vs y(s(x-u)) dx 


(2.49) 
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Let us denote the dilation of vCx) with a factor s by 
y (x) = Vs y (ax) 

A wavelet transfornt can be rewritten as inner products in 

L^R) 

Uf.(a.u) = < f(x), y^Cx-u) > (2.50) 

It thus corresponds to a decomposition of f(x) on the 
family of functions u«R^‘ functions y^(x) have the 

same type as y(x), but have a support s times smaller. If wavelet 
y(x) and signal f(x) have real values, in order to reconstruct 
f(x) from its wavelet transform, the Fourier Transform y(<>>) of 
y(x) must satisfy 

r I V'(")| 

Cy = ds <+« (2.51) 

The function y(x) can be interpreted as impulse response of 

A 

a bandpass filter. Let us denote y (x) = y (-x). We can rewrite 

8 8 

the wavelet transform at a point u and a scale s as a convolution 
product with y (x) 

Uf (s,u) = f * y„(u) . (2.52) 

8 

A wavelet transform can therefore be viewed as a filtering 

A. 

of f(x) with a bandpass filter whose impulse response is y (x) . 

8 

The fourier transform of y (x) is given by 

8 

y (<o) = — y (<•>/«) (2.53) 

* ys 
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In oppoaitlon to a window Fourier tranaform which haa a 
fixed reaolution in the apatial and frequency donain, the 
reaolution of a Uavelet Tranaforn variea with the acale paraineter 
a. Since v»(x) ia real, | yC**) | * | vC"")!- Let ua the centre 

A 

of the paesinft band of 

J j = 0 (2.54) 

o 

Let €f be the rms bandwidth around 6> 

o 

* J* I (2.55) 

o 

A 

It ia clear that centre of the paaaband of ia sut and 

s o 

that its r.n.a bandwidth ia sow. On a lofiarithmic scale, the rns 

bandwidth of y (<■>) is the sane for all a « R . Hence a wavelet 
a 

transforia deconpoaea the signal into a set of frequency bands 
having a constant size on a logarithm acale. 

2 

Let o-^ be the standard deviation of (yCx)] around zero. One 
can also show easily that the wavelet VgC* “ Uq) haa an energy 
concentrated around u^ within a standard deviation In the 

frequency domain, we saw that its energy is concentrated around 
a6>- within a standard deviation so-^ . In the phase-space, the 
reaolution cell of this wavelet ia therefore equal to [UQ-ir^/a), 
+ C<y /a)] X [S6)„ - so' , at*- + J . As opposed to a window 

0 U U 6> U a> 

Fourier transform, the shape of the reaolution cell variea with 
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the scale s- This is illustrated in Fig. 2.7. When the scale s is 
siaall, the resolution is coarse in the spatial dontain and fine in 
the frequency domain. If the scale s increases, the resolution 
increases in the spatial domain and decreases in the frequency 
domain (Fig. 2.7). 




Fie- 2.7 In the phase-epace, the ahape o£ a wavelet reaolutlon 
cell dependa upon the acale. Uhen the acale increaaea, 
the reaolutlon increaaea in the apatial domain and 
decreaaea in the frequency domain. The aurface of all 
the reaolutlon cell ia the aame. 




CHAPTER - 3 


RESULTS 


3.1 OF WAVELETS 


The 

wavelet functions w(x) 

have 

been 

constructed 

from 

orthogonal 

basis functions. Amongst 

the 

known 

orthogonal 

basis 

functions 

are Linear Spline, Haar 

basis 

, Quadratic Spline 

, and 


Cubic Spline. 


Ue start with a function A such that A are an orthonornal 

on 


basis for V^. Since ^ c span 


|4^(2.-n)| 


that 




• I 


^ (2x-n) 


there exists C such 
n 


(3.1) 


Define then 


(3.2) 


W(x) = ^ (-1)" ^(2x+n) 

Then , m,n « constitute an orthonornal basis of 

2 

wavelets for L (R). 

Now we take up orthogonal basis functions one by one. 

3.1.1 Haar Basis (31 

The orthogonal function is given by 


^(x) 


1 0 i X < 1 

0 Otherwise 
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is also exprsssed as 

= 4^C2x) + <fr(2x-l) 
Since <^(x) = Cq^(2x) + Cj^<^(2x-1) 
Therefore, 


Cq =1, = 1. 

By equation (3.2), 

V»(x) = ^(2x) - 0(2x-l) 


Hence 


1 

VCx) = -1 

0 


0 £ X £ 1 
1/2 :S X < 1 
Othervise. 


3.1.2! LlxMAr S^llnsr 131 


The function is given as 


X 

^(x) = 2-x 

0 


0 < X < 1 

1 < X i 2 
Otherwise. 


^(x) is also expressed as 

^(x) = j ^(2x) + ^C2x-1) + j ^(2x-2) 

1 1 
Hence ** * ^l""^ * ^2** 


(3.3) 


(3.4) 


(3.5) 
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By equation (3.2) 


vCx) = - Y <J^(2x+l) + ^(2x) - i ^(2x-l) 

Hence 


vCx) = -(X + i) 

. C3. - 
= Cj - X) 

= (^ - X) 

3. 1 . 3 (tuadratrlc Selina 


- ^ i X S 0 

os xi ‘ 
is xs 1 

1 S X s ‘ 


The quadratic apline function is 


given by 


<^(x) 


-2x^+6x-3 

(3-x)^ 


0 :£ X :£ 1 

1 :£ X :£ 2 

2 < X :£ 3 


= 0 


Otherwise. 


( 3 . 6 ) 


0(x) is also expressed as 

^(x) = j ^(2x) + j ^(2x-l) + I ^(2x+2) + i 4^(2 x- 3) 
Hence 



By equation (3.2) 

y»(x) = j ^(2x+2) - j ^(2x+l) + -y ^(2x) + ^ ^(2x-l) 



Hence 


vCx) 


- i (2x + 2)^ 

2 1 
-5x - 4x - -J 

2 1 
lOx - 4x - j 

-lOx^ + 16x - ~ 

5x - 14x + Y~ 

(4 - 2x)^ 


-1 aJ X a: - 
- i i X i 0 

0 a: X 

~ ^ X 1 

1 :£ X ^ '|- 

3 

X £ 2 . 


3.2! DESIGN OF FILTERS 


h(n) i« the Impulse response o£ discrete filter H. 
expressed as 

h(*i) ® ^ J (x/ 2 ). 4 ^ (x-n) dx. 

3. 3. 1 Hamr Basis 


For the haar basis function, h(n)’8 has been calculat 
tabulated. 


Tabla 3.1 



3.3.3 Linear Saline 


For the linear spline function, h(n)’a have been cal< 


and are as follows. 




47 


Table 3.2 


n 

h(n) 

0 

0.25 

1 

0.416 

2 

0.25 

3 

0.041 

-1 

0.041 


3.2.3 Ctuadrat,lc Spline 

For the quadratic apline function, h(n)'a have 
calculated and are aa follova : 


Table 3. 3 


n 

h(n) 

0 

0.637 

1 

-1.729 

2 

-3.479 

3 

-5.625 

4 

0.1208 

5 

0.004 

-1 

1.995 

-2 

0.004 


been 
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3. 2. <4 Cubic Saline 


The values o£ h(n)'s are as follows [2] : 


n 

h(n) 

n 

h(n) 

0 

0.542 

6 

0.012 

1 

0.307 

7 

-0.013 

2 

-0.035 

8 

0.006 

3 

-0.078 

9 

0.006 

4 

0.023 

10 

-0.003 

5 

-0.030 

11 

-0.002 


3.3 TRANSFER PlJ»K:TI(»iS - FILTERS 

The transfer function of the filters H & G corresponding to 
Haar basis and linear spline has been found out and are plotted 
(Figs. 3.1 and 3.2). The nature of the filters is highpass and 
lowpass respectively. 

3.4 i:HEC<»fP0SIT10N AND REC(M«STRllCTION BY WAVELET TRANSFfMtM 

5.^.4 

C13 Exponential : 

Ue have taken a signal, which is as follows : 
f(x) = 2 exp [ - ] + 2 exp [ - ] 

This signal was saapled and was then decomposed using the 
iterative method mentioned in Fig. 2.1. After the decomposition, 
the signal was reconstructed back by using reverse iterative 
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method mentioned in the Fig. 2.2. These procedures were carried 
out by taking the filters H & G corresponding to cubic spline and 
linear spline into account (Figs. 3.3 and 3.4). 

In both the cases, the reconstructed waveform follows the 
original waveform i.e. the reconstructed waveform has also an 
exponential decay. 

In case of cubic spline, the reconstructed waveform has some 
oscillations before they die down to zero. This may be due to the 
interpolation which is done during the reconstruction of the 
signal . 

In case of Linear Spline there is a sudden dip at the 
beginning in the reconstructed waveform. This may be due to the 
nature of the orthogonal basis function which decides the impulse 
response of the filters. 

C113 Sinusoidal s 

Ue have taken a sinusoidal signal which is as follows 
f(x) = 2 exp (-2ix) + exp (-ix) + 0.5 exp (-0.5 ix) 

The signal was sampled and was decomposed and 
reconstructed by utilising filters of cubic spline and linear 
spline (Figs. 3.5 and 3.6). In both the cases, although the 
reconstructed waveform does follow, the original waveform, but 
there are irregularities, the reasons again may be due to the 
interpolation. 
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3. 4. a DECOMPCKSITION AND RECONSTRUCTION OF SIGNALS WITH NOISE 
USING WAVELET TRANSFCWN 

Signals, both exponential as well as sinusoidal were 
analysed by adding white gaussian noise. The purpose o£ the 
adding noise and then reconstructing the signal from its 

decomposed components is to see whether or not there is an 

improvement in the signal quality. 

Cl> Exponential : 

Ca> Cubic Spline i The decomposition and reconstruction 

of noisy exponential signal was done by using filters 
based on cubic spline. It was observed that there was 
no improvement in the signal (Fig. 3.7). 

Cb3 Linear Spline t Uhen filters of linear spline were 
used, there was improvement of the signal quality. 
The reconstructed curve, although has a dip at the 
beginning, is smoother than the original noisy signal 
(Fig. 3.8). 

C113 SlnuBoldal a 

In case of a noisy sinusoidal signal, the signal 
quality is not improved compared to the original noisy signal 
using cubic spline as orthogonal basis function (Fig. 3.9 ). 
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3. 4r. 3 I%;C<MPOSITI<»f AND RECONSTRUCT! CM OF SI^AL WITH SPIKES 
C13 Exponential 

The exponential signal with very closely spaced spikes 
was taken for analysis. This signal was decomposed and 

reconstructed with cubic spline and linear spline. It was found 
out that in the reconstructed signal the spikes were resolved 
<Figa. 3.10, 3.11 and 3.12). 

Cii3 Sinusoidal 

The sinusoidal signal with closely spaced spikes was 
also taken for analysis. It was decomposed and reconstructed 
using cubic spline. The reconstructed signal showed the spikes to 
have been resolved (Fig. 3.13). 

3.S DECOMPOSITION AND RECONSTRUCTION BY LAPLACIAN PYRAMID SCHEME 

The above waveforms have also been decomposed and 

reconstructed by using Laplacian Pyramid scheme. 

Cl> Ebqpx»nent.lal : 

The exponential signal when decomposed and 
reconstructed with the help of Laplacian Pyramid Scheme, does 
follow the original curve. But some spikes of large amplitudes 
are seen inbetween. These may be due to the interpolation which 
is carried out while decomposing as well reconstructing the 
signal. Otherwise the curve Is smooth (Fig. 3.14). 
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C11I> Sinusoidal : 

This curve by undergoing decomposition and 
reconstruction through Laplacian Pyramid scheme, does follow the 
original curve. But again some spikes are visible in the 
reconstructed waveform. These may be due to the reasons as stated 
above for exponential curve (Fig. 3.15). 

3.6 ENERGY PLOTS OF IKTAIL SICMAL 


spikes 

The energy of the detail signal for the sinusoid with 


at various resolutions such 

as 0 

.5, 0.25, 0.125, 

and 

0.0625 

was 

found out and were plotted 

to 

observed that 

there 

were 

two 

distinct spikes occurring in 

all 

cases (Figs. 3 

.16. 

3.17, 

<0 

vH 


and 3.19). 

3.7 SEPARATION OF EXPfMENTIALS £5} 

The separation of exponentials is an age old problem. The 
idea is to find out the number of exponentials which make a set 
of data points. Also it is required to estimate the amplitude and 
time constant of the same exponentials. This problem has been 
dealt with at various stages. One of them is modelling approach. 
In this thesis, this problem has been taken up for tackling from 
the wavelet transformation point of view. The basic motivation 
for looking into the problem in this context is that when we 
carry out wavelet transform, the transformed data carry more 
information compared to the original data. 
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Uhen ve examine the decomposition and reconstruction of the 
exponential signal with the help of linear spline function we 
find that there is a dip at the beginning of the reconstructed 
signal before it converges to zero. That dip may be intepreted as 
an exponential with large time constant. Ue reconstructed the 
exponential with large time constant taking linear spline into 
account. Then we kept adding to it another exponential by varying 
its amplitude and time constant and these curves were decomposed 
and reconstructed. Now these were matched with the reconstructed 
curve of the original exponential function. Then we tried to 
find out which is the one closest to the original curve. 

The original exponential curve was 

fCx) = 2 exp 1^-0. 01 xj + 2 exp |^-0.1 xj 

The exponential which was kept fixed was 2exp(-0.01 x). The 
exponential which was varied was 2exp(-0.1 x). It was varied in 
the following manner : 

C13 AmplltrUde consbanb and bima consbanb varying i 

The amplitude was kept constant at 2, but the time 
constant were varied as 0.025, O.OS and 0.09 respectively 
(Figs. 3.20 and 3.21). 

Cil3 Amplitude varying and constant time constant s 

In this case the time constant was kept at 0.1, but the 
amplitude of the exponential was varied as 0.5, 1 and 1.5 
respectively (Figs. 3.22 and 3.23). 
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It waa observed that in the first case, the reconstructed 
exponential varied front each other and also from the original 
one. The exponential with time constant 0.09 was the closest to 
the original. But in the other case, when the amplitude is varied 
it was found out that all of them follow each other very closely 
and are very closed to original curve. 



CHAPTER - 4. 


CONCLUSIONS 


In this 

paper, we reviewed 

the 

application of 

wavelet 

decomposition 

and reconstruction 

of 

signals from 

several 


viewpoints. Ue covered exponential and sinusoidal signals showing 
that such a decomposition and reconstruction can be applied to 
such signals, lie described the mathematical properties of these 
decomposition. Ue first analysed different types of orthogonal 
basis functions and the corresponding wavelets. The properties of 
a window Fourier transform have been reviewed and a comparison 
has been drawn between the transform and the wavelet transform. 

Ue described the classical pyramidal mult iresolut ion 
algorithms and the wavelet approach to multiresolution 
decompositions. The difference of information between 
approximation of a function at two different resolutions is 
computed by decomposing the function into a wavelet orthonormal 
basis. Ue also explained the relationship between the orthonormal 
wavelet and quadrature mirror filter. The application of the 
wavelet decomposition and reconstruction to signals such as 
exponential and sinusoidal was carried out. It was found out that 
the linear spline is more applicable to exponential signals where 
there is a smooth reconstruction and also elimination of noise. 

As far as the separation of exponentials is concerned the 
problem was tackled from the wavelet point of view. It was found 
that while doing decomposition and reconstruction with linear 
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spline one o£ the exponentials which constitutes the original 
data is exposed. The other one is found out by trial and error 
method. This problem may be tackled by fitting a polynomial curve 
of exponentials to the reconstructed data. 
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